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Abstract 

Monitored Drift Tube (MDT) cliambers will constitute the large majority 
of precision detectors in the Muon Spectrometer of the ATLAS experiment 
at the Large Hadron Collider at CERN. For commissioning and calibration 
of MDT chambers, a Cosmic Ray Measurement Facility is in operation at 
Munich University. The objectives of this facility are to test the chambers 
and on-chamber electronics, to map the positions of the anode wires within 
the chambers with the precision needed for standalone muon momentum mea- 
surement in ATLAS, and to gain experience in the operation of the chambers 
and on-line calibration procedures. 

Until the start of muon chamber installation in ATLAS, 88 chambers built 
at the Max Planck Institute for Physics in Munich have to be commissioned 
and calibrated. With a data taking period of one day individual wire positions 
can be measured with an accuracy of 8.3 fim in the chamber plane and 27 /im 
in the direction perpendicular to that plane. 



1 Introduction 

The large majority of the precision detectors of the muon spectrometer of the AT- 
LAS experiment at the Large Hadron CoUider (LHC) wiU be Monitored Drift Tube 
(MDT) chambers. These chambers consist of 2 multilayers, each built of 3 or 4 
layers of densely packed drift tubes mounted on a support frame made of aluminum 
(compare figure^. The drift tubes — which are also made of aluminum — have 
an outer diameter of 3 cm and a wall thickness of 400 /xm; in the middle of each 
tube a gold-plated W-Re anode wire of 50 /xm diameter is stretched. The anode 
wire is positioned only at the two tube ends by precision end-plugs; the sag of each 
chamber can be adjusted to follow the sag of the wire along the tube. The drift 
tubes are operated at a pressure of 3 bar absolute with an Ar:C02=93:7 gas mixture 
at a gas gain of 2 x 10*. The area covered by a complete chamber ranges from 1 m^ 
to 11 m2. A total of 1194 MDT chambers will be used in ATLAS [T]. 



Figure 1: A schematic view of an ATLAS MDT chamber. Also indicated are the 
optical paths of the geometry monitoring system (in-plane alignment). 

In collaboration with the Max Planck Institute (MPI) for Physics, Munich Uni- 
versity (LMU) is responsible for the construction of 88 MDT chambers of type BOS 
(Barrel Outer Small — the second largest type used in the ATLAS barrel region) 
which consists of 432 drift tubes, arranged in 2 x 3 layers of 72 tubes, and has a size^ 
of 3.9 m X 2.2 m x 0.5 m. The chambers are built with high mechanical precision at 
MPI |2] and are then commissioned and calibrated at the Cosmic Ray Measurement 
Facility of Munich University. 

The main objectives of the Cosmic Ray Facility are to test the chambers, to 
map the positions of the tube layers and of individual anode wires within the cham- 
bers, and to gain experience in the operation of the chambers and the calibration 
procedures. 

An important design goal of the ATLAS muon spectrometer is the standalone 
measurement capability of muon momenta up to 1 TeV with a relative error of less 
than 10%. This leads to the requirement that the anode wire positions within the 
MDT chambers must be known to within 20 /xm Q . 

2 The Cosmic Ray Measurement Facility Setup 

A schematic view of the Cosmic Ray Measurement Facility is shown in figure |21 
The whole setup is located in a climatized hall to ensure stable environmental 

^Some chambers will be narrower or have cut-outs — realized with shorter drift tubes — to 
accommodate support structures of the ATLAS detector or fcedthroughs to inner components. 
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Figure 2: A schematic view of the Cosmic Ray Measurement Facihty. The individ- 
ual components are described in the text. 



conditions during the measurements. We define a right handed coordinate system, 
the X coordinate in the direction along the wires, y perpendicular to the wires in 
the chamber plane, and z perpendicular to the chamber plane pointing upwards. 

Three MDT chambers are operated simultaneously in the facility. Two of them, 
the so-called reference chambers, encompass the third chamber — which is to be 
tested — and provide track information for cosmic ray muons. The anode wire 
positions of the reference chambers have been mapped with a precision of 2 /.im 
with an X-ray tomograph at CERN 

A coincidence of the three scintillator layers in the hodoscopes above and below 
the MDT chambers provides the trigger. The upper hodoscopc consists of one 
layer of scintillator counters, the lower of two staggered layers displaced by half the 
width of a counter. Only muons with a momentum greater than 600 MeV pass the 
34 cm iron absorber above the lower hodoscopc and can generate a trigger signal. 
The trigger logic is divided into 5 segments along the x coordinate to limit the 
inclination of the muons in the x-z plane. The arrangement of the scintillator 
bars allows the measurement of muon tracks in the x direction with a precision of 
8 cm. The time of the muon transition is determined by the lower hodoscopc with 
a precision of better than 800 ps. The average trigger rate is about 70 Hz. 

The streamer tubes at the bottom of the setup are used together with the lower 
reference chamber to determine the multiple scattering angle in the y-z plane in the 
iron absorber. This angle is used in combination with the angle between the track 
segments in the two reference chambers to estimate the momentum of the muon. 

Two contact free alignment systems are used to continuously monitor the po- 
sitions of all MDT chambers: The position of the reference chambers relative to 
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each other is measured with an optical ahgnment system consisting of eight RasNik 
sensors @] , the position of the test chamber with respect to the upper MDT cham- 
ber with eight capacitive sensors |S] . Both systems have a precision of better than 
5/im 

In addition to the chambcr-to-chamber ahgnment, the internal chamber geom- 
etry is also monitored with integrated RasNik sensors (compare figure^. 



3 Drift Time Spectra 

The position measurement in the MDT chambers is based on the measurement of the 
drift time of the electrons (from ionization by the incident muon along its trajectory) 
to the anode wire. The drift time is determined with respect to the external trigger 
system and translated to the impact radius (distance from the anode wire) via the 
so-called r-t relation (see section 0J); the drift tubes are only sensitive to the track 
coordinates in a plane perpendicular to the anode wire, not along it. The drift 
of the electrons depends on the gas composition, the gas density, and the electric 
field applied (i.e. the operating voltage), and is therefore characteristic of the gas 
mixture and the operating conditions. 

After correcting for the signal propagation times along the anode wire and in the 
electronics and cables, the distributions of drift times are used as a first check that 
each tube is functioning as expected and that the detector response is homogeneous. 
For the simple case of uniform illumination of a tube, the density of hits at a given 
drift time is proportional to the drift velocity at the corresponding radius. 

Events with more than one charged particle passing through the setup would 
lead to ambiguities in the hodoscope time measurement. Therefore, exactly one hit 
in the upper hodoscope and one pair of hits in overlapping counters in the lower 
hodoscope are demanded, which rejects 15% of all triggers (see table ^1. 

The edges of the drift time spectra are parameterized [J] with the functions 

F (t) = Pa -\ _t for the rising edge and (1) 

1 + e"V 

G (t) = p„i + '^-t for the trailing edge of the spectrum, (2) 

1 + 

as shown in figureOl The parameters po and p„i correspond to the rates of accidental 
hits, Aq to the height of the spectrum, and to to the mid-time and Tq to the steepness 
of the rising edge of the spectrum, which is determined by the resolution of the drift 
tube near the anode wire. Correspondingly, is the middle of the trailing edge 
of the spectrum and Tm its steepness, which is determined by the decrease in pulse 
height near the tube walP. The parameterization with a„it + Am accounts for the 
slope of the spectrum before its trailing edge that is due to the decrease of the drift 
velocity with increasing radius. 

The maximum drift time, defined as r = (i,„ + T„i) — to, is characteristic of the 
drift properties of a tube. The distribution of the maximum drift time is shown 
in figure 2Ji for the tubes of one MDT chamber. As part of the chamber quality 
control programme, it is checked that all tubes have the same drift time spectrum 
characteristics within errors, compare figure ^ja. 

^The small pulse height for muons passing a tube close to the tube wall is due to their short 
path inside the tube, which leads to smaller primary ionization. 
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normalized drift time spectrum 




t/ns 



Figure 3: A typical drift time spectrum, normalized to one entry. The fitted 
functions arc explained in the text. 




Figure 4: a) Distribution of the maximum drift time r = (tm + Tm) — to for the 
tubes of a tested chamber, b) The corresponding pull distribution. 



4 The Relation between Drift Time and Drift Ra- 
dius 

The relation between the measured drift time and the corresponding drift radius {r-t 
relation) is determined without the help of external spacial detectors and without 
assumptions on uniform tube illumination. After the check described in section |21 
it is assumed that all tubes of an MDT chamber have the same r-t relation. 

Events with two hits at small drift times (|i — to| < 15 ns) in one MDT chamber 
are selected. By assuming that the muon went through the two corresponding 
anode wires and that the trajectory of the muon is a straight line, the impact radii 
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in the other tubes of the chamber can be calculated. This procedure yields drift 
times for a discrete set of drift radii. The number of such points — typically 16 — 
is determined by the geometry of the chamber and the angular acceptance of the 
trigger. In addition to these points, the drift times corresponding to r = and to 
the inner tube radius are determined by the fit to the drift time spectra described in 
sectional The r-t relation is then approximated by a linear interpolation between 
the measured values. 



5 Track Reconstruction 

Within each MDT chamber a muon track segment is reconstructed from the mea- 
sured hits. A track segment is assumed to be a straight line in the y — z plane'^: 

y{z) = m„z + bn, (3) 

where the index n specifies the chamber. The parameters to„ and 6„ are determined 
for each chamber by minimizing the function 

racket '^drift,i) / .\ 

where the index i runs over all tubes contributing to the track segment, rtrack,i is 
the distance of the track segment to the wire, rdrift,i is the measured drift radius, 
and Ui is the resolution of the tube for the given radius. At a test beam measurement 
using an external reference tracking system of silicon strip detectors the intrinsic 
single tube resolution was determined to be around 100 /im |S]. 

The tubes contributing to a track segment are selected by choosing one tube 
with a hit in each layer of the chamber, which is a sample of six tubes. From the 
impact radii in the two outermost tubes of this sample four possible track segment 
candidates are derived. If the difference between the distance of such a candidate 
to the wire and the measured drift radius does not exceed 2 mm in any of the other 
tubes, a track segment is fitted using this sample. If more than one track segment 
is found, the one with the smallest is used. If no track segment is found, the 
procedure is repeated with subsamples of five tubes. 



6 Chamber Alignment 

The MDT chambers in the Cosmic Ray Facility are aligned mechanically with a 
precision of about 100 /im. A much more precise estimate of their relative posi- 
tions is then obtained by comparing the muon track segments reconstructed in each 
chamber. The optical and capacitive alignment systems are used to verify that the 
setup does not move significantly during the data taking period. In the future, the 
alignment information will be included in the analysis to correct for time-dependent 
position shifts of the chambers. 

The positions of the two reference chambers are measured relative to the test 
chamber. 

• A shift Vy of a reference chamber with respect to the test chamber results in a 
systematic shift between the track segments in the reference and test chamber 
and can be obtained from the average over all events, 

Vy = {brcl - &tost) , (5) 

^The third coordinate (x) is only used to correct for signal propagation times and for the wire 
sag and is taken from the hodoscopc information. 
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where ^lef and &tcst sltc the track segment parameters defined in equation |j3Jl- 
They arc determined with a precision of 1 fim. 

• A shift Vz resuhs in a relative shift of the track segments depending on the 
slope m. It is determined from a linear fit to the distribution of 6iof — fetcst 
versus m with an error of 10 /im. 

• A shift Vx along the x coordinate can be ignored, as the MDT chambers are 
not sensitive to such a shift along the direction of their anode wires. 

• The tilt angle a of a reference chamber around the x axis is given by the 
systematic deviation of the two track segment slopes as 

a = {rUrel - mtest) ■ (6) 

A precision of 10^® is achieved. 

• To determine the angles P and 7 corresponding to rotations of the reference 
chamber around the y and z axes, the x coordinate measurement of the ho- 
doscope is used. The shifts Vy and Vz are calculated for three one meter wide 
sections in x. From the dependences Vz{x) and Vy{x), the angles /3 and 7 are 
obtained with precisions of 8 x 10~® and 10~^, respectively. 

In the further analysis, the parameters of track segments in the reference cham- 
bers are corrected according to the measured reference chamber positions. 



7 Energy Estimation 



The energy of cosmic muons traversing the Cosmic Ray Measurement Facility is 
estimated from multiple scattering — which depends on the muon momentum — 
in the chambers and the iron absorber. 

If TOref.u designates the slope of the track segment reconstructed in the upper 
reference chamber and miof,i the slope of the track segment reconstructed by the 
lower reference chamber, the distribution of Am = m^cf, u — Wref.i will be the wider 
the lower the muon energy. The standard deviation aAm{E^) has been determined 
as a function of the muon energy by means of a Monte Carlo simulation and can be 
parameterized as [3] 



<^Am{E^j) ~ (Jo 



/eOOMeV 
V E^ 



(7) 



with fToo = (2.0 ± 0.2) • 10-4, = (1.001 ± 0.008) • 10"^ and a = 1.04 ± 0.01. 

Similarly, the deviation Aj/ of the position measured with the streamer tubes 
from the prediction obtained from the track segment in the lower reference chamber 
depends on the muon energy: For its width 'J\y{E^), a Monte Carlo simulation 
yields the values CToo = (4.4 ± 0.1) mm, ctq = (183 ± 2) mm, and a = 1.44 ± 0.01. 

The probability density / for measuring deviations Am and Ay is then given by 



f{E, 



1 



27r CTAm{Efj) (JAy{Efi) 



cxp 



Am 



(JAm{Ef^) 



Ay 



(JAyiE^ 



(8) 



For a given track, the value of E^ that maximizes f{E^) is an estimate of the muon 
energy. This estimator is biased and has a limited resolution but is still useful for 
a selection of high momentum muons. 



8 Measurement of Wire Positions 

In this section, the method for the measurement of anode wire positions is presented. 
To determine wire positions, the muon track segments in the reference chambers 
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are extrapolated into the test chamber, and the drift radii measured in the test 
chamber are compared with the track predictions. 

The Cosmic Ray Measurement Facility is particularly sensitive to displacements 
6y of the wires in the chamber plane. From this, a shift of entire tube layers in the 
y direction and the mean spacing g between the wires within a layer can also be 
determined. 

Wire displacements 5z in the plane perpendicular to the chamber are accessi- 
ble via tracks with different inclinations m (cf. equation (jSJ)- This measurement 
is less precise than the 6y determination owing to the limited angular acceptance 
of the trigger and the angular distribution of cosmic muons. During chamber pro- 
duction, z displacement and tilts of entire tube layers are more difficult to control 
than individual wire displacements Sz with respect to the layer. Therefore, the z 
displacement for each tube layer is also measured. 

The measurements are performed at both ends of the chamber by selecting 
muons which passed the tested chamber within one meter of the corresponding 
end of the chamber using the hodoscopc information. Because the wires are only 
supported at the tube ends and since the sag of the anode wires is known from 
wire tension measurements these two measurements completely determine the wire 
position. 

8.1 Event Selection 

Events with one reconstructed charged particle are selected as described in sec- 
tion For the further analysis, the event selection proceeds in two steps: First, 
the reconstructed track segments are subjected to a set of cuts, and second, the 
individual hits in the tubes must pass certain quality criteria. 

The presence of a reconstructed track segment with at least 5 hits is required 
in each of the three MDT chambers. The track segments must also roughly match, 
i.e., |6ref,u — bre[,i\ < 4mm and |mref,u — "^ref,i| < 15mrad. This ensures that the 
track segments are not affected by tubes in which delta rays''' distorted the drift time 
measurement. The remaining events are called good tracks. For the determination 
of wire positions, tracks with high muon momentum are desirable to limit the effects 
from multiple scattering. Thus, the estimated muon energy is required to be larger 
than 2.5 GeV. The efficiencies of these track cuts are given in tabled 

Tubes in the test chamber with hits generated by delta rays would affect the 
measurement of the wire position, and are therefore rejected: It is demanded that 
the tube have contributed to a track segment in the test chamber. Also, a loose cut 
on the difference between the radius prediction from the reconstructed track seg- 
ments and the measured drift radius is made. For the track segment reconstructed 
in the test chamber this cut is at 0.7 mm, while it is at 1.5 mm for the predicted 
radius from the reference track segments extrapolated into the test chamber^. The 
efficiencies of these cuts are given in table El They are very close to 100% and thus 
do not bias the wire position measurement. 

8.2 Determination of the Wire Positions 

The reconstructed reference chamber track segments are extrapolated into the test 
chamber and the drift radius measured in the test chamber is compared with the 

* Energetic primary electrons knocked from the tube walls or gas atoms by the incident charged 
particle. Delta rays can lead to smaller measured drift times if they are emitted in the direction 
of the anode wire. 

^The larger value for the reference chamber track segments is duo to the long extrapolation 
distance into the test chamber. 
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Figure 5: A sketch explaining the method for the measurement of wire displace- 
ments. A muon (dashed line) passing a tube with a wire that is displaced from 
the nominal position by 6y in the y direction and Sz in the z direction will lead 
to a measured drift radius rdrift, while rrcf is predicted by the reference chambers 
relative to the nominal wire position. 



prediction from those track segments. The weighted average r,.ct of the drift radius 
predictions from the two track segments 

_ (l/gg)rref,u + (l/gf)r,cfj 

is taken, where each uncertainty ai (i = u, 1 for the track segment in the upper and 
lower reference chamber) depends on the extrapolation distance from the reference 
chamber and the estimated muon momentum (to account for the uncertainty from 
multiple scattering). The value of ai is about 200 j-im for the reference chamber 
closer to the tube hit and ranges from 350 to 800 ijm for the farther chamber. 

If a given wire is displaced by 6y in the y direction and 6z in the z direction, the 
measured drift radius rdrift in the tube will show a systematic deviation Ar from 
the predicted reference radius, 

Ar = Tdrift - r^of = Sy - mSz + 0{m'^) , (10) 

as shown in figure [S] 



8.3 Determination of the Wire Position Perpendicular to the 
Chamber Plane 

The displacement 6z of an individual anode wire is determined from the slope in a 
linear fit to the distribution of Ar vs. m. 

The z displacement of an entire tube layer and its tilt corresponding to a rotation 
about the x axis are determined from the distribution of measured Sz values of all 
tubes as a function of their y coordinate^. A linear fit to this Sz distribution is 
performed separately for each layer. The offset and slope of the fitted line are taken 
to be the z displacement and tilt of the layer, respectively. 

®Here, y is taken to be the nominal wire position. Because of the large chamber width, this 
approximation does not bias the result. 
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In the further analysis, each wire is assigned the measured z position of the layer 
at the nominal y wire position. The single wire position measurement is used for 
wires where it deviates from the layer position by more than 3 standard deviations. 

8.4 Determination of the Wire Position in the Chamber Plane 

Using the z position measurement described above, the displacement by of an anode 
wire is obtained as 

(Ar') = (r;,f - rdrift) ^ by ~ {m)5z' « 5y , (11) 

where r^^j denotes the predicted radius corrected for the wire z position and the 
residual displacement 5z' after correction is of the order of the resolution of the 
z position measurement. The average of the track slopes (|(?7i)aii hits |) all tubes is^ 
0.06, the largest value of |(m)aiihits| for a single tube being 0.15. Therefore, the 
systematic error introduced in the last approximation of equation 111() is typically 
much smaller than 15% times the resolution on the z position, which does not limit 
the 5y measurement. 

In the chamber plane, the nominal wire positions can be described by a grid of 
the form 

y{n) =yo + ng , (12) 

where n denotes the index of the tube in the layer, yo the offset of the layer position 
relative to the other layers in the same chamber, and g the mean distance between 
two neighbouring wires. Fitting the function in equation H12|) to the distribution of 
the wire positions of one tube layer yields the offset yo and the mean distance g for 
that layer. 

8.5 Test of the Method 

In order to test the method described above, one of the few chambers which have 
been scanned by the X-ray tomograph at CERN has been used as a test chamber. 
Two tomograph scans of the chamber were performed at a distance of 30 cm from 
each chamber end. Therefore the positions of the anode wires® are known with a 
precision of 2 ^m. 

The analysis uses 28 hours of data taken with this chamber. During this time 
1.7 million events per meter along the anode wires {x direction) have been recorded. 
The number of events after each step of the event selection and the number of tube 
hits in the test chamber used for the wire position measurements are given in tables^ 
andm 

All measurements at the Cosmic Ray Facility are in very good agreement with 
the tomograph results. 

In figure El the 6z measurements for single wires performed at the Cosmic Ray 
Facility are compared with the results of the tomograph scans. The distributions of 
the difference between the cosmic ray measurement and the tomograph scan have 
widths of 23 /xm and 27 /im. These widths are dominated by the resolution of the 
Cosmic Ray Measurement Facility. 

A better precision can be obtained for the z displacement of entire tube layers. 
A comparison of this measurement and the values for the tilt angle of the tube layers 
around the x axis (the parameters which are more difficult to control during the 
MDT chamber assembly) with the tomograph results is shown in figures and |H1 

^The deviation from is caused by a slight asymmetry of the trigger efficiency in the y direction 
during data taking. 

*Only typically 80% of the wires are visible in the tomograph scans because of the support 
structure of the chamber. 
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selection cut 


number of 
events after cut 


efficiency 
of cut 


all events 


1.65 • 10*^ 




hodoscopc hits from one particle 


1.40 • 10*^ 


85% 


good tracks 


0.84 • 10^ 


60% 


estimated energy > 2.5 GcV 


0.34 • lO^ 


40% 



Table 1: Number of events per meter along the wire after each cut in the event 
selection (second column) and the percentage of events which pass the cut (right 
column) . 



selection cut 


number of 
tube hits after cut 


efficiency 
of cut 


all tube hits 
hits on tracks 
residuum cuts 


2.04- 10^ 
1.90- 10^ 
1.87- 10^ 


96% 
98% 



Table 2: Number of tube hits in the test chamber per meter along the wire (as 
used for the wire position measurement) after each cut in the event selection (second 
column) and the percentage of hits which pass the cut (right column) . 



Side A Side B 




Figure 6: Comparison of the Sz measurement for single wires at the Cosmic Ray 
Facility with the tomograph scans. Note the different scales in the two upper plots. 
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Figure 7: Measurements of 8z for entire tube layers. The triangles denote the 
cosmic ray measurement, the dots the X-ray tomograph data. 
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Figure 8: Measurements of the tilt of layers. The triangles denote the cosmic ray 
measurement, the dots the X-ray tomograph data. 
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Figure 9: Comparison of the Sy measurement for single wires at the Cosmic Ray 
Facihty with the tomograph scans. 



Here, we achieve a precision of 4.4 /im for the layer shift and 17 x 10^^ for the tilt 
angle. 

The comparison of the Sy measurement with cosmic rays with the tomograph 
results shows that the difference between the two measurements has an RMS of 
8.5 fiui at cither end of the chamber (see figure^. The precision of the Cosmic Ray 
Measurement Facility is therefore 8.3 /im (after subtraction of the uncertainty on 
the X-ray tomograph measurement). 

The measurements of the layer offset yo and the wire grid constant g are com- 
pared in figures ITUl and ITTI A precision of 1.8 /im for the layer offset and 0.15 jum 
for the grid constant is achieved. 



9 Conclusions 

In the Cosmic Ray Measurement Facility, ATLAS MDT muon chambers are tested 
and calibrated. The Cosmic Ray Measurement Facility is capable of measuring the 
positions of the individual anode wires of a chamber with a precision of 8.3 /im in 
the chamber plane and 27 /im in the direction perpendicular to that plane after 28 
hours of data taking. 

In addition to the wire positions, the offsets of the tube layers relative to each 
other in the chamber plane and their grid constants have been determined with 
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Figure 10: Measurements of the layer shift in y. The triangles denote the cosmic 
ray measurement, the dots the X-ray tomograph data. 
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Figure 11: Measurements of the deviation of the grid from the nominal grid. The 
triangles denote the cosmic ray measurement, the dots the X-ray tomograph data. 
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precisions of 1.8 /xm and 0.15 /^m, respectively. The positions of the tube layers 
in the direction perpendicular to the chamber plane and the angles for a rotation 
around an axis parallel to the wires have been measured with precisions of 4.4 fim 
and 17 x 10^^, respectively. 

At a rate of up to two tested chambers per week, all 88 BOS chambers which 
are produced in Munich will be tested until the planned start of muon chamber 
installation in ATLAS. 
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